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We study dynamical pliase transitions at temperature T = in a fully frustrated square Josephson 
^\. , junction array subject to a driving current density which has nonzero components ix, iy parallel to 

Qs ■ both axes of the lattice. Our numerical results show clear evidence for three dynamical phases: a 

pinned vortex lattice characterized by zero time-averaged voltages {vx)t and {vy)t, a "plastic" phase 
in which both {vx}t and {vy)t are nonzero, and a moving lattice phase in which only one of the time- 
average voltage components is nonzero. The last of these has a finite transverse critical current: if a 
current is applied in the x direction, a nonzero transverse current density iy is required before {vy)t 
becomes nonzero. The voltage traces in the moving lattice phase are periodic in time. By contrast, 
the voltages in the plastic phase have continuous power spectra which are weakly dependent on 
frequency. This phase diagram is found numerically to be qualitatively unchanged by the presence 
S ' of weak disorder. We also describe two simple analytical models which recover some, but not all, 

O ' the characteristics of the three dynamical phases, and of the phase diagram calculated numerically. 
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I. INTRODUCTION 



There has recently been considerable interest in dynamical phases in superconductors. Such interest has been 
'"O ■ stimulated by much evidence (both theoretical and experimental) that vortex lattices, as well as other ordered systems, 
C can exhibit various types of transitions from one type of dynamical phase to another, as a function of controllable 

external parameters, such as driving force and temperature. In the case of Type-II superconductors in a magnetic 
field, the driving force which acts on the vortex lattice is the applied current. 

In disordered superconducting films and crystals, there now appear to exist at least three distinct phases, as a 
function of driving current and temperature K-M- At low driving current, the vortex lattice is pinned, and typically 
exhibits a glass-like order because of the random pinning centers which prevent the lattice from moving. At intermedi- 
ate driving currents, the glass phase is depinned, and starts to move; this motion is thought to occur inhomogeneously 
f— ^ ' in disordered superconductors (that is, some of the vortices move through random channels, while other vortices re- 
\^ , main pinned). In this phase, the vortex system is usually said to exhibit plastic flow. Finally, at high driving currents, 
f^ ■ an event resembling a phase transition occurs, and the vortex system reverts back to an ordered state which closely 
0^ , resembles a moving crystalline phase. Most workers now believe that this third phase, especially in two dimensions, 
^^ ■ lacks long-range crystalline order but instead exhibits hexatic or smectic order, modified by the static disorder of the 
system [n,B[. This high-current phase is also thought to exhibit a finite transverse critical current, i. e., once moving 
parallel to the driving force, the lattice requires a finite driving force perpendicular to the average direction of motion 
Ch ' before acquiring a nonzero transverse velocity. The details of all these phases, as well as the transitions between them, 
r^ , remain the subjects of much experimental and theoretical investigation. 
r^ ' Recently, there has been evidence that some similar behavior is to be found in ordered systems, both in two and 

Q , three dimensions (2D and 3D). Extensive numerical studies have been carried out by Nori and collaborators, who 
O • have found a very complex structure of ordered and disordered (sometimes called "plastic" ) phases as a function of 
^ [ the relative density of vortices and pinning sites, and of the direction of applied current relative to the axes of the 
periodic pinning lattice [p|jlO[ . Some of the observed structure can be described in terms of nomenclature frequently 
used for incommensurate dynamical systems, such as devil's staircases, Ar'nold tongues, etc. These workers use a 
^ , standard model whereby the vortices interact with each other via a pairwise potential involving a modified Bessel 
■ ■ ' function, and with the pins by a suitable short-range pinning potential. A similar model had been used previously 
to study vortex lattices driven through a disordered array of pins uM- This earlier calculation revealed a somewhat 
simpler structure with three phases apparent: a pinned vortex lattice, a moving plastic phase, and a moving vortex 
lattice with a finite transverse critical current. 

In this paper, we investigate the possibility of dynamical phase transitions in a widely investigated model system: 
the overdamped square Josephson junction array in a transverse magnetic field. This study is complementary to 
the work mentioned above, in that the phase angles of the complex superconducting order parameter are explicitly 
included as degrees of freedom. Although such ordered Josephson junction arrays have, of course, been extensively 
studied |12[, little attention has been paid to the possibility of plastic phases in these materials. Strongly disordered 



arrays have recently been shown numerically to have both plastic and moving-lattice phases separated by a phase 
boundary which appears to exhibit critical phenomena associated with a diverging correlation length |13|| . Weakly 
disordered driven periodic systems are also predicted to have a nearly periodic temporal order at large values of the 
driving parameter |14| . 

In order to search for possible dynamical phases, we study this system in a regime in which a current is applied with 
nonzero components parallel to both the x and y axes of the array. We consider only the so-called fully- frustrated case 
with exactly 1/2 flux quantum per plaquette of the square array. Even in this ordered array, we find three phases. 
These are (A) the pinned vortex lattice, (B) a moving lattice, and (C) a plastic phase occurring at currents between 
A and B. Phases A and B have been much discussed previously, but phase C seems not to have been observed in an 
ordered lattice. The plastic phase C is characterized by aperiodic voltage noise and an absence of long range vortex 
order. By contrast, the moving lattice B is characterized by voltage signals which are periodic in time, and, most 
strikingly, by a finite transverse critical current. 

In an attempt to understand our numerical results, we also present some simple analytical models. One model is 
based on the assumption that the dynamics can be characterized in terms of periodically repeating 2x2 unit cells 
driven by a current applied at an angle to the cell axes. A second model considers the motion of a single vortex through 
the periodic potential formed by the lattice of Josephson-coupled grains. Both models give rise to distinct regimes 
in which none, one, or both of the voltage components parallel to the cell axes are nonzero. Thus, in particular, a 
nonzero transverse critical current is present in both analytical models. 

The remainder of this paper is organized as follows. In Section II, we briefly review our calculational model. Our 
numerical results are presented in Section III. Section IV is an interpretive discussion which includes the two analytical 
models mentioned above. Our conclusions are described in Section V. 

II. MODEL 

In all our simulations, we take as the starting point the well-known dynamical equations for an overdamped 
Josephson junction array in a transverse magnetic field. We assume that each Josephson junction is resistively 
shunted, and we neglect the loop inductance; that is, we assume that the magnetic field generated by the currents is 
negligible compared to the applied magnetic field. The equations then take the form 

lab = Ic,ab sin {(j)a - <j)b - Aab) + TTS— ZH (<^a ^ A ^ ^ab) , (1) 

leHab at 






lab = IT- (2) 



Here Ic,ab is the critical current of the junction connecting grains a and 5, Rab is the shunt resistance of that junction, 
and /|^* is the external current fed into the ath grain. We assume a constant, uniform external field B — Bz 
perpendicular to the array, and make the gauge choice A = Bxy. The phase factor 

Aab=^ [ A-ds (3) 

5'o J 

is then easily expressed in terms of / = Ba^/^o, the frustration or flux per plaquette of dimension a x a, where 
^o — hc/2e is the flux quantum. 

The equations of motion can be put in dimensionless form with the definitions iab = lab/Ici ic;ab = Ic,ab/Ic, and 
gab = R/Rab, and the natural time unit t = ^J^j , where Ic and R are a typical critical current and a typical shunt 
resistance. The result of this substitution is 

iab = ic,ab sin (^a ^ 4>b - Aab) + ^dab-T: (0a " (t>b) , (4) 



T.'-'> = ^r■ (5) 

b 

Combining these equations yields a set of coupled differential equations which is easily reduced to matrix form and 
solved numerically pq ]. In our work, we employed a fourth-order embedded Runge-Kutta integration with variable 
time step [na. 



In most of our calculations, we have considered a lattice without disorder: all shunt resistances Rab = R (or gab — 1), 
and all critical currents ic,ab = 1 on a square array. Our array is driven by nonzero current densities in both the x 
and y directions: a current i^ is fed into each grain along the left-hand edge of the array, and extracted from each 
grain on the right-hand edge, and a current iy is similarly injected into each grain on the bottom edge and extracted 
from each grain on the top edge of the array (see Fig. nl). 

Starting from randomized initial phases, we integrate these equations of motion over an "equilibration" interval of 
100t-5000t, followed by an averaging period of 100t-2000t. Typically, the external currents i are ramped up or down 
(in steps of 0.001 to 0.1) without rerandomizing the phases. We calculate the spatially averaged but time-dependent 
voltage difference w(i') = V{t')/NRIc (where t' = t/r is the dimensionless time) between the input and output edges 
of the array in both the x and y directions, as well as its time-average (v)t, in both directions. In some regions of 
the phase diagram, these voltages appear to be periodic in time, as revealed by an analysis of the power spectrum of 
the voltage. In other regions, as described below, this power spectrum reveals that the voltage is aperiodic. In some 
of our simulations, we also tracked the number and motion of vortices in the array. The vortex number in a given 
plaquette a is an integer n„_Q defined by the relation 



n. 



D,a 



-^^(0a-0fc-A,b)-O,±l, (6) 



where the sum is taken clockwise around the a*'' plaquette, and each phase difference is restricted to the range [— tt, tt]. 
Our calculations are carried out exclusively for a frustration / = 1/2, i. e., an applied magnetic field equal to one 
flux quantum for every two plaquettes. In a square Josephson array, the ground state of this field is the well-known 
checkerboard vortex pattern, shown schematically in Fig. |l[ Our simulations reproduce this pattern. 

III. RESULTS 

The central results of our calculations arc summed up concisely in Fig. 0, which shows the "dynamical phase 
diagram" for two square Josephson junction arrays with no disorder at / = 1/2, driven by two orthogonal currents 
ix and iy. We find three different phases: a pinned vortex lattice (time-averaged voltages Vx = Vy = Q), & plastic flow 
regime {vx = 0, Vy> Q or Vy = 0, Vx > 0), and a moving vortex lattice {vx > 0, Vy > 0). The calculated boundaries 
are shown in Fig. g. 

The results of Fig. g were obtained by two different methods. In the flrst procedure, the longitudinal current ix 
was flrst ramped up from zero to a flnite value, with iy held at zero. Next, the transverse current iy was ramped up 
at flxed ix- To determine the phase boundaries, we simply searched for the currents at which the time averages {vx)t 
or {vy)t (or both) became nonzero. In the second method, we ramped up ix and iy simultaneously, holding the ratio 
iy/ix flxed. Both methods generally gave similar phase boundaries. Likewise, we found little indication of substantial 
hysteresis, except in determining the boundary between phases A and C. In this case, if the integration time is too 
long, the system tended to jump abruptly back and forth between the two phases. 

We refer to the phase in which both voltages are nonzero as a "plastic" phase, by analogy with a similar phase 
exhibited by vortices in systems with quenched disorder. In this phase, v{t') is apparently non-periodic in time; 
the corresponding voltage power spectrum is only weakly frequency-dependent (see below). Although we use the 
nomenclature of plastic, we have not checked that the vortex motion in this phase is inhomogeneous (i. e. that only 
some vortices are in motion while others remain pinned in this phase). Such inhomogeneous motion is thought to 
occur in disordered systems pT| . By contrast, in the driven lattice phase, where only {vx)t or {vy)t is nonzero, the 
power spectra of the voltage traces are sharply peaked at a fundamental frequency and its harmonics. Further, while 
one voltage is always nonzero, the other voltage is periodic, only averaging to zero over a cycle. We interpret this 
behavior as representing a vortex lattice being driven transverse to the larger of the two current components. 

Various minor numerical difficulties sometimes interfered with the calculations, but they could usually be overcome. 
For example, spurious voltage jumps were occasionally observed during these calculations; these jumps (unlike genuine 
jumps) could generally be eliminated by changing the initial conditions, the integration time, or the direction of current 
ramping. Only those voltage jumps which appeared at the same position on the phase diagram in different runs were 
deemed to be genuine. From the occurrence of these jumps, however, we conclude that at certain points on the phase 
diagram there are a several metastable dynamical states which have similar energies. The occurrence of such states 
may suggest a first-order transition across the phase boundary, at least in the ordered system. 

In Fig. pi we show time-dependent voltage traces at several points in the plastic and moving-lattice phases. The 
traces are plainly very different in the two phases. In the moving-lattice phase, the voltage traces are evidently 
periodic in time. By contrast, in the plastic phase, the voltages, both in x and y directions, are obviously aperiodic. 
Another striking feature is apparent in the moving lattice phase. In this phase, as noted above, there is a non-zero 



time averaged voltage only along one of the two directions, even though current is applied along both the x and y 
directions. Despite this feature, there is a finite time- dependent voltage in the y direction, which averages to zero, 
and which is periodic like Vx{t'). 

To learn more about the harmonic content of these voltage traces, we have also calculated the voltage power spectrum 
P{ujt), using the non- normalized Lomb method for variable-time-step data p7|]: 
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where to is defined by 



X;^ cos2 uj {t.j - to) J2j sin^ w (tj ~ to) 



y^ . sin 2ujti 

tan(2^io) = ^^ — i. (8) 

2^ cos2iLjtj 



Here, lo is the angular frequency, the ij's are the times at which the voltage is recorded, Vj = v{tj), and v is the 
arithmetic average of the Wj's. We have carried out this calculation for several points in the plastic and moving-lattice 
phases as indicated by triangles in Fig. 2(b). The results are shown in Fig. for both Vx and Vy. Clearly, the 
noise confirms that the voltage in the moving-lattice phase is periodic in time, while that in the plastic phase has a 
continuous spectrum which is relatively weakly dependent on frequency. 

The moving-lattice phase is characterized by a fundamental angular frequency loq. It is readily shown that loq is 
related to the time-averaged voltage drop across the array, {V)t, by 

coo = 2e{V)t/h (9) 

This relation is consistent with the widely-accepted egg-carton picture of the vortex lattice in the moving phase [Q . 
In this picture, the vortex lattice is viewed as a collection of eggs moving in a potential similar to an egg-carton, 
consisting of a periodic distribution of wells on a square lattice (each well lying at the center of a plaquette formed by 
four grains). During one period, the vortex lattice moves by one row through the egg-carton potential. Since there is 
a phase slip of 27r between two opposite edges of the array each time a vortex crosses the line joining those two edges, 
one can readily deduce the above relationship. 

Fig. 0(a) shows that a finger of the driven lattice phase is interposed between the pinned and plastic flow phases in 
a 10 X 10 array. This finger appears to be a finite-size artifact, because it is absent in the phase diagram for a 20 x 20 
array shown in Fig. ||(b). We believe that the phase diagram of Fig. g(b) is likely to persist in an A^ x A^ lattice even 
at very large N. Thus, the ordered array at / = 1/2 has three phases: pinned lattice, moving plastic phase, and 
moving lattice. 

Since a realistic Josephson lattice is certain to have some disorder, we have also carried out a limited number 
of calculations for an array at / = 1/2 with weak disorder in the critical currents. Specifically, we assume that 
the critical currents are independent random variables uniformly distributed between 0.9/c and 1.1/c. The resulting 
phase diagram, for a single realization of disorder, is shown in Fig. 0. It is calculated using the same techniques 
as for ordered arrays. Once again, we see clear evidence of three phases: pinned lattice, plastic phase, and moving 
lattice. These have characteristics similar to those in the ordered case. For example, the moving lattice phase has a 
finite transverse critical current ic±, which goes to zero near the phase boundary. In this 10 x 10 sample, there is an 
even larger finger of B phase interposed between A and C than there is in the ordered array; once again, we assume 
that this finger disappears in larger arrays. For strongly disordered square arrays at several different field strengths, 
a phase diagram resembling ours, though without the interpolated finger, has been found by Dominguez |13|| . We 
have also calculated the voltage power spectra for the phases B and C; they resemble those of Fig. 4 in that P{ujt) 
has peaks at multiples of a fundamental frequency in the moving lattice, while the spectra in the plastic phase are 
continuous and only weakly dependent on frequency. 

IV. SIMPLIFIED ANALYTICAL MODELS 

In a further effort to understand the behavior found numerically, we have considered two simple analytical models. 
In this section, we give a brief description of the models used. As will be seen, while each can reproduce some of the 
numerical results, neither generates all the details of the simulations. 



A. Four-Plaquette Unit Cell 

Our first analytical model is a slight generalization of an approach previously used by Rzchowski et al p9| | . They 
consider the dynamics of a fully frustrated N x N array of overdamped resistively-shunted junctions. To treat this 
system analytically, they assume that the dynamical state is simply a periodic repetition of a square four-plaquette 
unit cell. It has long been known H^] that the ground state of the fully frustrated lattice has such a unit cell, 
corresponding to the checkerboard vortex pattern shown in Fig. |l|. If we maintain this assumption of periodicity, the 
equations of Ref. [|l9[ are readily extended to to the case of currents applied at an angle to the plaquette edges. The 
resulting equations take the form 

j + a + fi + d = TT{mod2TT) (10) 

d(3 dS d'y da 

_sm/3---smJ--+sm7+^+sma+- = (11) 

d'y dd 

_+sin7-^-sina = /,„,, (12) 

dp d5 ^ 

— + sm /? - — - sm 6 = Itot.y (13) 

where a, /3, 5, and 7 are the four inequivalent gauge-invariant phase differences describing the bonds of the 2x2 
primitive cell (cf. Fig. 1 of Ref. fl^), t' = t/r is a dimensionless time, and Itot,x and Itot,y are the total bias currents 
in the x and y directions per 2x2 superlattice cell (in units of the single-junction critical current). 

We have solved these equations numerically, first reducing the system to three variables and then employing the 
same integration algorithm described above. In comparing with the equations for our first set of simulations, note 
that the quantity Itot.a = 2ia{a = x,y) is twice the current injected into each boundary grain. The resulting phase 
diagram is shown in Fig. |6[ As in the previous phase diagrams calculated in this paper, there are regions (denoted A, 
B, and C) in which none, one, or both of the time-averaged voltages Vx and Vy are nonzero. Despite the reduction 
in number of variables and the enforced symmetry of this simulation, the voltages in regime C are still aperiodic in 
time, with continuous power spectra at most points in regime C. A representative power spectrum for Vx{t') is shown 
in the inset to Fig. g. It is calculated for a point in region C indicated by a triangle. 

When the current is applied along the x axis, the critical current is close to the analytically computed p^ value 
ix = \/2 — 1- By contrast, when ix » 1, the boundaries of region B asymptotically approach the line iy = 1. In 
general, we conclude that this simplified version of the dynamics has some, but not all, features of the large array. In 
particular it does have a region of "plastic flow" at large ix and iy. However, it fails to reproduce the region of plastic 
flow interposed between phases A and B and seen in our larger-scale simulations. 

B. Single Vortex in an "Egg-Carton" Potential 

Our second analytical model is even simpler. It refers, not to an entire lattice of vortices, but to a single vortex 
moving in the "egg-carton" potential formed by the lattice. According to Lobb, Abraham, and Tinkham |1§|| , a single 
vortex can be viewed, to a good approximation, as moving in a potential of the form 



Vix,y) = -Vo 



2ttx\ /27ru\ 
cos I 1 -I- cos I 1 



(14) 



where a is the lattice constant of the array and Vq is the depth of the potential felt by a single vortex. This potential 
has minima at r = {x,y) = (nio, 7120) where n = 0,±1,±2, .... [The grains of the Josephson lattice, in these 
coordinates, are located at ((ni + i)a, (71-2 + |)a), and correspond to maxima of the vortex potential.] 

The current-voltage characteristics of this model are readily calculated. The Magnus force on a vortex due to an 
external current density J may be written (taking z as the direction perpendicular to the array) 

Fext = ^oz X J/c, (15) 

where J, in this two-dimensional system, represents a current per unit length. In the steady state, this force has to 
be balanced by two other forces: the gradient of the egg-carton potential energy and the frictional force experienced 
by the vortex moving through the lattice. This condition may be written 

Fext - VV{y) - ryr = 0, (16) 



or explicitly, in component form, 

fix = Vosm (17) 

c a \ a J 

■qy = Vo sm . (18) 

c a \ a / 

Given the vortex velocities, the electric fields may be written down by using the relation between vortex velocity 
and electric field. Specifically, the voltage drop between any two points Pi and P2 is given by 

AV^i2 = 2TTnv^L — , (19) 

2e 

where vj_ is the component of vortex velocity perpendicular to the line joining Pi and P2, ^ is the distance between 
Pi and P2, and n is the vortex number density per unit area. Here we have used the fact that the phase difference 
between Pi and P2 changes by 27r every time a vortex crosses that line, and have also availed ourselves of the Josephson 
relation between voltage and phase. We therefore deduce the following expressions for the components of electric field: 

E. = ^ (20) 

hni 
E, . - — . (21) 



Eqs. ( |l7|) and (18) are identical in form to the equations for single Josephson junctions, with x/a and y/a playing 
the role of the Josephson phase, 27rVo/a the role of the critical current, and — $0 Jj,/c and ^qJj./c the roles of the 
driving currents. In view of eqs. (O) and dlq), we deduce that the time-averaged voltage drop in the i direction 
{i — x,y) becomes nonzero when |$o>^i/c| > 27rVo/a. As has been shown by ref. [1^], Vq is related to the critical 
current /j, of an individual Josephson junction by 

Vo ~ 0.22^. (22) 

Collecting all this information, we obtain the simple phase diagram shown in Fig. M for the voltage drops arising from 
motion of a single vortex in a square array. Once again, there are three regimes, denoted A, B, and C, where none, 
one, or both of the time-averaged voltage drops (wx)* and {vy)t are nonzero. However, the phase diagram is simplified 
by the fact the {vx)t and {vy)t are independent of one another, depending only on i^ and iy respectively. The phase 
boundaries correspond to the vertical and horizontal lines Vx ~ 0.11 and Vy ~ 0.11. 

The results of Fig. [7| are obviously oversimplified compared to the numerical diagrams of Fig. 2. In addition to 
the other obvious differences, the voltages in region C are not chaotic, but are instead just the superposition of two 
independent voltages in the x and y directions, each of which has its own fundamental frequency and harmonics of 
that fundamental. The two fundamentals may, of course, be incommensurate depending on the values of {v,j;)t and 
{vy}t. 

V. DISCUSSION AND CONCLUSIONS 

The present calculations clearly show that a fully frustrated array of overdamped Josephson junctions exhibits at 
least three dynamical phases as a function of the two orthogonal driving currents ix and iy: a pinned vortex lattice, 
a plastic phase characterized by a continuous power spectrum for both v^ and Vy, and a moving vortex lattice with 
only one of the voltages Vx and Vy nonzero. In this last phase, the power spectrum, at least for the limited lattice 
sizes we have investigated, contains only harmonic multiples of a fundamental frequency. Weak disorder in the critical 
currents appears not to change this phase diagram greatly. 

Regarding our phase diagram, it is natural to ask whether the boundaries between the different phases are analogous 
to first-order phase transitions. Although we have no conclusive evidence, our numerical results suggest that they 
may indeed be first-order, rather than continuous, at least for the ordered lattices. In support of this conjecture, we 
note the occasional occurrence of hysteresis in our simulations, and of discontinuous jumps between one phase and 
another near the phase boundaries. There is also little evidence that any quantities, such as the strength of the voltage 
noise, diverge near the phase boundaries, as might be expected of a continuous phase transition. In this respect, these 
transitions differ somewhat from those seen in strongly disordered lattices. 



For these ordered arrays in which the vortex lattice is commensurate with the underlying Josephson array, one might 
have expected other types of commensurate-incommensurate transitions as the angle between the applied current and 
array symmetry axis is varied. Such "magic angle" effects, with a multitude of commensurate and disordered phases, 
are seen in other models in which vortex lattices are driven through periodic pinning arrays ||9|,|l0[|. In the present 
case, we have seen no clear evidence of any phases other than the three shown in our phase diagram. Possibly, such 
additional phases would appear if we studied larger arrays in which such delicate effects would be more stable. On 
the other hand, our open boundary conditions may constitute such a strong perturbation on the periodic array that 
such commensurability effects would be suppressed even for very large lattices. 

Finally, we comment briefly on the power spectra found in our simulations. In the moving lattice phase B, the 
voltage power spectra appear to contain only multiples of a fundamental frequency. Such a power spectrum represents 
an array which is phase-locked, and hence would radiate power only at multiples of that fundamental frequency. In 
this case, phase-locking clearly occurs between rows as well as along rows of junctions parallel to the voltage. Such 
phase locking must be present because it is required in order to produce the observed relation between {vx)t and the 
fundamental frequency. Furthermore, such locking survives weak disorder in the critical currents. The effectiveness 
of a magnetic field in producing phase locking in square arrays has been discussed previously i2l| , |22|| . The present 
work provides additional evidence that a field / = 1/2 is effective in producing phase locking, even when the applied 
current has nonzero components parallel to both array axes |23|| . 

In conclusion, we have numerically investigated the dynamic phases of fully-frustrated square Josephson junction 
arrays driven by two independent, orthogonal currents. We find three phases: stationary lattice, driven lattice, and 
plastic flow. These phases appear also in weakly disordered arrays. We find that our numerical results can be partly 
understood by two simple analytical models, which show some of the features of the full simulation. 
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FIG. 1. Schematic of the geometry for calculating IV characteristics. Currents i^ = I^/Ic and iy = ly/Ic are injected into 
each grain (shaded circles) on left-hand and lower edges of the square array, and extracted from right-hand and upper edges, 
as shown. Time-dependent but spatially averaged voltage drops Vx{t') and Vy{t')a.Te calculated across the array in the x and y 
directions. Crosses denote vortex locations, as defined in the text, for the ground-state checkerboard pattern at / = 1/2. 

FIG. 2. Calculated zero-temperature "dynamical phase diagram" for (a) 10 x 10 and (b) 20 x 20 square Josephson arrays 
with currents applied in both the x and y directions. Region A: pinned vortex lattice {{vx)t ~ {vy)t ~ 0). Region B: moving 
vortex lattice phase {{vx)t ~ 0, {vy)t > 0, or {vy)t = 0, {vj:)t > 0). Region C: moving plastic vortex phase {{vj:)t > 0, {vy)t > 0). 
The phase boundaries were determined numerically at the calculated points, as described in the text; the smooth lines are 
interpolations between these points. The phase diagram is assumed symmetric about a 45° line; diamonds denote points which 
were calculated for ix > iy; circles, for iy > ix. The sliver of phase C between A and B in (a) is a finite-size artifact, as shown 
in (b) where it is seen to be absent in a 20 x 20 array. Open triangles in (b) indicate currents {ix, iy) at which the voltage 
traces and power spectra of Figs. H and W were calculated. 

FIG. 3. Time-dependent voltage traces Vx{t') and Vy{t') (upper and lower parts of each panel) for several points in the phase 
diagram of a 20 x 20 array ai f — 1/2. (a) [ix = 0.351, iy — 0.079] and (b) [ix ~ 405, iy = 0.294], represent points in the plastic 
phase C; (c) [ix = 0.508, iy = 0.210] and (d) [ix = 0.469, iy = 0.287] are points in the moving-lattice phase B. As described in 
the text, Vx{t') represents the difference between the average voltages on the left-hand and right-hand edge of the array; Vy{t') 
is the difference between the average voltages on the top and bottom edges of the array, i' = t/r is a dimensionless time; r is 
a natural time unit defined in the text. 

FIG. 4. Calculated voltage power spectra PioJT) for the voltages Vx{t') and Vy{t') of Fig. H (a) - (d). Here co is the angular 
frequency. In each panel, the power spectrum for Vx lies above that for Vy. 

FIG. 5. Calculated phase diagram for a 10 x 10 array at / = 1/2 with one realization of weak disorder. Rather than all 
being equal, the critical currents are chosen at random from a distribution which is uniform in the interval (0.9/c, 1.1/c). 

FIG. 6. Calculated phase diagram for a square array at / = 1/2, obtained by solvingthe overdamped dynamical equations 
for a 2 x 2 unit cell which is assumed to be repeated throughout the lattice [eqs. hol-llS]. The phase diagram is shown as a 
function of the driving currents ix and iy in the x and y directions. As in previous phase diagrams, there are regions where 
none, one, or both of the two time-averaged voltages {vx)t and {vy)t are nonzero. Inset shows power spectrum of Vx{t') for the 
point indicated by the open triangle {ix — 1.348, iy = 1.115). 

FIG. 7. Calculated phase diagram for a single vortex moving in the egg-carton potential produced by a square array of 
overdamped Josephson junctions. The notation is the same as in Fig. H. Once again, there are regions where none, one, or 
both of the time averaged voltages {vx)t and {vy)t are nonzero. 
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